Spin Polarized and Valley Helical Edge Modes in Graphene Nanoribbons 



O 

(N 



Zhenhua Qiao, 1 Yugui Yao, 2, 1 Shengyuan A. Yang, 1 Bin Wang, 3 and Qian Niu 1 

1 Department of Physics, The University of Texas, Austin, Texas 78712, USA 
'^Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
3 Department of Physics, The University of Hong Kong, Hong Kong, China 

Inspired by recent progress in fabricating precisely zigzag-edged graphene nanoribbons and the 
observation of edge magnetism, we find that spin polarized edge modes with well-defined valley index 
can exist in a bulk energy gap opened by a staggered sublattice potential such as that provided by 
a hexagonal Boron-Nitride substrate. Our result is obtained by both tight-binding model and first 
principles calculations. These edge modes are helical with respect to the valley degree of freedom, 
and are robust against scattering, as long as the disorder potential is smooth over atomic scale, 
resulting from the protection of the large momentum separation of the valleys. 
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Introduction. — The appearance of edge states is one 
of the most peculiar phenomena in solid state systems. 
They are often connected to topologically non-trivial 
bulk properties, e.g. non-zero Chern numbers in quan- 
tum Hall systems [l|, Q, or odd Z2 numbers in time- 
reversal invariant topological insulators [3j. The edge 
states in quantum Hall systems are robust against all 
kinds of disorders and interactions [J], while those in 
the latter systems can survive scatterings that preserve 
the time reversal symmetry 0, The edge states in 
graphene with zigzag terminations belong to a different 
category. Such states connect the two different valleys 
K and K' projected along the edge direction and their 
presence is dictated by the bulk topological charge [fj. 
It is of great interest to utilize these unusual states for 
various applications @, [1| ■ 

Recently, zigzag-edged graphene nanoribbons have 
been fabricated with precision by unzipping carbon nan- 
otubes [ii-HH . Without electron-electron interaction, the 
edge states form a completely flat edge band connecting 
the two valleys with large momentum separation 12 , 13[ . 
When interaction is taken into account, due to the sin- 
gular density of states, spins on the edge become spon- 
taneously polarized resulting in an edge ferromagnetism 
0, fl3 | , w hich has been confirmed by a recent experi- 
ment [151] • The spin polarized edge states are disper- 
sive in momentum space, making them useful for current 
transport. Unfortunately, without a bulk gap, the effect 
of edge states would be overwhelmed by the contribution 
from the bulk states. 

In this Letter, we show that for zigzag-edged graphene 
nanoribbons spin polarized dispersive edge states can ex- 
ist and remain robust in a bulk gap opened by a staggered 
sublattice potential. Such potential can be realized, for 
example, by a hexagonal Boron-Nitride substrate. The 
edge modes are then tied to the conduction or valance 
band edges, and with spontaneous spin polarization due 
to interaction effects, one spin branch of the edge modes 
is pushed into the bulk gap forming a conducting chan- 
nel in the bulk insulating state. We show that such spin- 



polarized states within the bulk gap remain robust under 
scattering with correlation length longer than the lattice 
constant. We also perform first principles calculations to 
further support our predictions. 

Tight-Binding Model. — Figure [I] (a) illustrates the 
schematic setup of a zigzag-edged graphene nanoribbon 
in the presence of a staggered sublattice potential. A 
tight-binding Hamiltonian that incorporates phenomeno- 
logically the edge spin-polarization can be written as: 



H 
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where c|(cj) is the electron creation (annihilation) oper- 
ator on site i, and a z is the z component of Pauli ma- 
trices. The first term is the nearest neighbor hopping 
with t being the hopping energy. The second term repre- 
sents the effect of edge ferromagnetism which stems from 
electron-electron interaction, and Mq is a phenomeno- 
logical parameter put in by hand at present stage, and 
its value will be determined later from the first-principles 
method. The last term corresponds to the staggered A/B 
sublattice potentials: f7j=A/2 for sublattice A (o), and 
Ui=— A/2 for sublattice B (•). In our analysis of tight- 
binding model, we measure the energy e, magnetization 
Mo, potential A, and disorder strength W in units of the 
hopping energy t. 

The energy spectrum of graphene nanoribbons with 
zigzag termination can be numerically obtained by diag- 
onalizing the Hamiltonian H(k) in momentum space for 
each crystal momentum k along the edge direction which 
we take as x-axis. Figure [2] shows the evolution of the 
energy spectrum as functions of staggered sublattice po- 
tential A/2 and edge magnetization Mo for fixed ribbon 
width N = 800 (about 852 A). For clarity, only the edge 
states from the upper boundary are shown. Panel (a) 
shows the doubly-degenerate flat-bands connecting the 
two Dirac points K and K' . We observe that a bulk 
energy gap A = 0.4 is opened by inversion symmetry 
breaking due to the staggered sublattice potentials. 
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FIG. 1: (a) Diagram of zigzag-edged graphene nanoribbons. 
A (o) and B (•) sublattices are subjected to staggered po- 
tentials. Magnetization Mo only exists at the edge atoms in 
the grey regime. TV labels the ribbon width, (b) Schematic 
diagram of edge states propagation in a two-terminal mea- 
surement geometry for Fermi energy in the bulk gap of 
Figs. [2Ib)-(d). (c) First principles calculation of the spin- 
polarization P for the carbon atoms of the supercell of zigzag 
graphene on top of the hexagonal Boron Nitride substrate. 
P= (n t -nj.)/(n t 



Figures [5] (b)-(d) plot the band structures when Mq is 
taken to be 0.6, 1.0, and 1.4, respectively. We find that, 
due to the different degrees of localization of the states in 
edge band [f| , the magnitude of the energy splitting of the 
edge bands is fc-dependent: the spin- up edge band bends 
upward, while the spin-down edge band bends downward. 
This makes edge band dispersive hence capable of con- 
ducting charge current. Moreover, for a fixed sublatticc 
potential A/2=0.2, we find that along with the increasing 
of the edge magnetization Mo from 0.6 [see panel (b)] to 
1.0 [see panel (c)], the spin-down edge band gradually ap- 
proaches the bulk valence band, and eventually touches 
and emerges into the bulk valence band (at Mqc±.1A). 
This creates gapless edge modes tied to each valley, which 
is similar to the finding in Ref. @ except that the edge 
modes here are spin-polarized. 

From the energy spectra, we observe that the two edge 
states in the bulk energy gap propagate along opposite di- 
rections. Due to the same flat-band origin, they are local- 
ized at the same upper boundary. We also note that the 
edge states in the gap have well-defined valley index K 
or K' . The situation is schematically shown in Fig.QJb): 
the edge states with same spin but different valley indices 
propagate oppositely along the same boundary. This can 
be natually termed as a spin-polarized quantum-valley 
Hall state. 

Robustness of Spin- Polarized Edge Mode. — From the 
above analysis, we notice that for the fixed bulk gap, the 
spin polarized edge state is gapped for weak edge mag- 
netization, while the edge state becomes gapless when 
the edge magnetization approaches a critical Mq. These 
edge states in the gap provide conducting channels for 
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FIG. 2: Evolution of the band structures of zigzag edged 
graphene nanoribbons at fixed width TV = 800 (only exhibit- 
ing the edge states from the upper boundary), (a) When the 
staggered sublattice potential A = 0.4 is applied, an bulk gap 
is opened, and the flat-bands are doubly-degenerate; (b)-(d) 
The edge magnetism is further included with Mo = 0.6, 1.0, 
and 1.4, respectively. The flat bands become spin-split: spin- 
up edge band bends upwards, while spin-down edge band 
bends downwards. 



spin-polarized transport. However, to be useful for prac- 
tical applications, they need to be robust against impu- 
rity scattering. In the following, we shall investigate the 
robustness of the edge state in the presence of impuri- 
ties, and show that both gapped and gapless edge states 
are robust against scattering due to the large momentum 
separation between the valleys K and K' . 

It is known that the impurity scattering in graphene 
mainly comes from the long-range Coulomb scatter- 
ers [16[. We assume that the impurity potential Vi at 



each site i takes a Gaussian form [171 ]: 



Vi = ^ iv j x exp(- 



2xf 



(2) 



where the summation is over all sites, Wj is the local 
disorder strength at site j, and £ is the correlation length. 
We define an effective disorder strength W from £ and 
w [H: 



W = w x (£ 2 + 1) 



(3) 



The numerical simulations are performed within the same 
setup of Ref. [20j including only the left and right semi- 
infinite leads. The two-terminal conductance can be cal- 
culated from the Landauer-Buttiker formula 



21|: 



G = -Tr[r R GT L G a ] 



(4) 



where G r ' a are the retarded and advanced Green's func- 
tions of the central disordered region. The quantities 
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FIG. 3: Averaged conductance (G) (in units of e 2 /h) versus 
effective disorder strength W at fixed A=0.4, and Mo=0.6 for 
different Fermi energies e = — O.f (square), O(circle), O.f (trian- 
gle), respectively, (a) Long ranger disorder (5=4); (b) Short 
ranger disorder (5=0). 20 000 ensembles are collected for each 
set of parameters. 



Tl/t? are the line- width functions describing the coupling 
between the left/right lead and the scattering region, and 



can be obtained from Tp = i(E^ 



Here. 



2jp are 



the retarded/advanced self-energies of the P — th semi- 
infinite lead, and can be numerically evaluated using the 



22|. 



recursive transfer matrix method 

Figure [3] plots the sample-averaged two-terminal con- 
ductance (G) as a function of the effective disor- 
der strength W for three different Fermi energies e=- 
0.1, 0, 0.1, respectively. The edge magnetization is set 
to be Mo=0.6. Each data point represents average over 
20,000 sample configurations. Panel (a) is for the long 
range disorder case. We observe that, for all the three 
Fermi energies, the averaged conductances (G) are ro- 
bust against weak disorders, i.e. when FF<0.8, all the 
conductances are exactly quantized to be one in units 
of e 2 /h without conductance fluctuation. When W>0.8, 
we find that (G) of e=-0.1 is quickly destroyed first, and 
that of e=0.1 is the most robust one. This can be ex- 
plained from the band structure shown in Fig.^b). One 
can see that the two edge states for a fixed Fermi energy 
have a large momentum separation when the Fermi en- 
ergy is near the upper band bottom (e.g. e=0.1). The 
separation decreases when the Fermi energy is approach- 
ing the valence band top (e.g. £=-0.1). The large mo- 
mentum separation (on the scale of valley separation) 
suppresses the long range impurity scattering which al- 
lows only small momentum transfer. Panel (b) shows the 
averaged conductance as a function of the short range 
nonmagnetic disorders, with other parameters being the 
same as that in panel (a). We find that the edge states 
are very sensitive to the disorders and therefore easily de- 
stroyed by small disorder strengths. We also performed 
calculations for magnetic disorders (not shown here) and 




FIG. 4: (Color online) (a) Atomic structure of hydrogen- 
terminated zigzag-edged graphene nanoribbons on top of sin- 
gle layer hexagonal Boron-Nitride. The red square regime 
represents a supercell. Upper: side view; lower: top view, (b) 
Band structure for spin anti-parallel configurations between 
two zigzag edges. In the bulk energy gap A around 78meV 
(the narrow energy window of the projected bulk band struc- 
ture shown in light grey), only spin-up states exist, and a triv- 
ial gap is opened due to the electron-electron interaction with 
the same spin between the two zigzag edge boundaries, (c) 
Band structure for spin parallel case. Spin-up and spin-down 
states coexist in the gap. (d) When a weak voltage bias 0.27V 
is applied transversely, the upper (lower) edge states are up- 
wards (downwards) lift, leaving only the spin-down states in 
the gap. Bands in Red and green represent spin-up and spin 
down edge bands, respectively. 



the results are similar. Therefore, we conclude that our 
valley associated spin-polarized edge modes are robust 
against smooth disorder scattering. 

First Principles Calculations. — So far, we have in- 
vestigated the spin-polarized edge modes in a graphene 
nanoribbon model Hamiltonian with a staggered AB 
sublattice potential, and an edge-specific spin polar- 
ization put in by hand. In the following, from first 
principles calculations, we provide a concrete system 
for the realization of our models by placing a zigzag- 
edged graphene nanoribbon on top of a hexagonal Boron- 
Nitride substrate. In our calculations, we set the lat- 
tice constant to be a — 2.45A, and inter-layer distance 
d = 3.22A [11. Fi gure|U(a) illustrates the schematic 
configuration of the system. Here, we use Ni (N2) to la- 
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bel the width of graphene (Boron nitride), and Ni<N2- 
The single-layer graphene and Boron nitride are AB 
stacked with Nitrogen atoms on top of the hollow po- 
sition. All the outer-most atoms are saturated with Hy- 
drogen atoms, and we use the experimental values of the 
bond lengths: 1.17A (B-H), l.OlA (N-H), and 1.09A (C- 
H). The self-consistent ground state calculations were 
performed within the non-equilibrium Green's function 
coupled with the density- functional theory scheme [23| . 
and the local density approximation exchange-correlation 
potential (LDA-PZ81) was used(24|. 



In our calculations, we set N\ — 96, and N2 = 112. 
Panels (b) and (c) plot the band structures with spin anti- 
parallel/parallel configurations at the two zigzag bound- 
aries. The grey region shows bulk band structure region 
(when both N\ and N2 approach infinity): a bulk gap 
A around 78 meV is opened, which is slightly larger 
than 53 meV using VASP package [25] . In panel (b), 
we find that only the spin-up polarized edge states lie 
inside gap, which resembles the band structure of the 
tight-binding model. A small splitting 5 in the figure is 
due to the interaction between the states with the same 
spin on the two boundaries, and will decrease along with 
the increasing system width. In panel (c), we observe 
that both the spin-up and spin-down states coexist in 
the gap. From both panels (b) and (c), one can obtain 
that the magnetization of each outermost carbon atom 
is about Mo=0.287 eV. As shown in panel (d), one can 
apply an external transverse bias V = 0.27 V across the 
ribbons 7] to separate the mixed states and leave only 
the spin-down edge state in the bulk gap, which provides 
an efficient way to manipulate the spin-polarized edge 
states. 

Conclusion. — We investigate on the edge modes of 
zigzag-edged graphene nanoribbons in the presence of a 
staggered sublattice potential. We find that the edge 
states form spin-polarized conducting channels which are 
robust against smooth impurity potentials. Using first 
principles calculation methods, we provide a specific sys- 
tem which exhibits such spin-polarized edge modes by 
placing the zigzag-edged graphene nanoribbons on top 
of a hexagonal Boron-Nitride substrate. The realization 
of this valley associated spin-polarized edge modes will 
enable the application of the graphene-based spintronics 
and valley tronics devices. 
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